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I .  INTRODUCTION 


In  1914,  Dougall  treated  the  problem  of  equilibrium  of  an  elastic 
rod  with  the  use  of  a  general  solution  of  the  equations  of  elasticity. 
Dougall ?s  solution  is  expressed  in  terms  of  three  arbitrary  harmonic 
potentials.  In  1960,  Muki^  treated  elastic  equilibrium  problems  with  the 
use  of  a  general  solution  expressed  in  terms  of  another  set  of  potentials. 
Muki?s  potentials  comprise  an  arbitrary  biharmonic  function  and  an 
arbitrary  harmonic  function.  Our  purpose  here  is  to  establish  the 
relations  between  the  Dougall  potentials  and  the  Muki  potentials. 


II.  THE  EQUATIONS  OF  ELASTIC  EQUILIBRIUM-NOTATION 


Our  analysis  will  be  in  terms  of  cartesian  tensors.  (For  a  further 
discussion,  see  Appendix  A.)  We  deal  with  a  cartesian  coordinate  system 
x-p  X2,  Xg.  We  shall  adopt  the  convention  that  lower  case  Latin  indices 
will  take  the  values  1,  2,  3.  A  comma  followed  by  an  index,  say  i, 
denotes  partial  differentiation  with  respect  to  the  coordinate  corres¬ 
ponding  to  that  index,  say  x^. 

Dougall  writes  the  equations  of  elastic  equilibrium  in  the  form 


y  V  ui  +  (X  +  y)  A, i  =  0, 


(1) 


2  . 


where  u^  is  the  displacement  vector,  V  is  the  laplacian  operator;  i.e 


V  u.  =  u.  ,  .  . 
i  i  JJ 


(2) 


A  denotes  the  volumetric  expansion;  i.e.. 


A 


u 


i’i 


and  A  and  y  are  the  Lame  constants. 


(3) 


Muki,  on  the  other  hand,  wrote  Eqs.  (1)  in  the  form 


V 


+ 


1 

1  -  2v 


0 


(4) 


J.  Dougall ,  nAn  Analytical  Theory  of  the  Equilibrium  of  an  Isotropic 
Elastic  Rod  of  Circular  Cross - s ectionn ,  Trans .  Roy.  Soc .  Edinburgh , 
XLIX,  Part  IV,  (No.  17),  pp  895-978  (1914). 


2 

R.  Muki,  nAsymmetric  Problems  of  the  Theory  of  Elasticity  for  a  Semi- 
Infinite  Solid  and  a  Thick  Platen,  Progress  in  Solid  Mechanics,  I, 
pp  401-439  (1960). 
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where  v  is  the  Poisson  ratio,  and  Eqs.  (2)  and  (3)  hold.  In  order  to 
compare  their  solutions,  we  should  express  them  both  in  terms  of  the  same 
material  constants.  We  find  it  notationally  convenient  to  write  the 
equations  of  equilibrium  given  by  Eqs.  (1)  and  (4)  as 


V 


2 

u. 

i 


0 


where 


X  +  v 

2  (X  +  2y) 


(5) 


(6) 


In  the  forms  of  solution  given  both  by  Dougall  and  by  Muki,  one  of 
the  axes  plays  a  distinguished  role.  We  shall  take  this  to  be  the  X3  axis. 
This  is  done  because,  for  problems  involving  cylinders,  it  is  convenient 
to  take  a  coordinate  axis,  usually  called  the  z-axis,  to  be  in  the  direction 
of  the  cylinder  axis.  Thus  we  shall  identify  x-^ ,  X2>  and  X3,  respectively, 
with  the  x,  y  and  z  axis  of  a  cartesian  coordinate  system.  We  shall  then 
adopt  the  convention  that  lower  case  Greek  indices  will  take  the  values 
1,  2  and  we  shall  identify  the  index  3  with  the  index  z.  Thus, 


x3  =  Z 


u_  =  u 
3  z 


u. , _  =  u. ,  , 

1  3  1  z 


etc.  We  then  write  our  equations  in  expanded  form.  For  Eqs.  (5)  we 
write 


V2u 


2A 


a  2A  -  1  ’a 


A.  =  0  , 


(7) 


n2  2A  a  n 

v  uz  -  2a~-  r  4,z  =  0 


(8) 
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For  Eqs.  (2)  we  write 


V  u  =  u  ,_0  +  u  ,  , 

a  a  33  a  zz 


V  u  =  u  >nn  +  u  , 
z  z'00  Z^zz 


For  Eq.  (3)  we  write 


A  =  u  ,  +  u  , 

a  a  z  z 


(9) 


We  are  now  ready  to  introduce  the  potentials. 

III.  THE  DOUGALL  POTENTIALS 

^Dougall  introduces  three  potentials*  which  we  shall  denote  as  <f>,  0 
and  \p.  These  potentials  are  arbitrary  harmonic  functions,  i.e., 

v2<f>  =  o,  v2e  =  o,  v2^  =  o. 

Dougall  states  his  general  solution  as  an  arbitrary  linear  combina¬ 
tion  of  three  basic  solutions,  which  are  as  follows:  one  solution  is 
given  by 


u 


X  <(>,  ,  u  = 

Y  zz  z 


-  X  (f>, 

y  yz 


2(X  +  2y)  ~ 

A  +  y  9’z 


(10) 


which  are  his  Eqs.  (2).  A  second  basic  solution  is  given  by 


Ui  =  9’i 


which  is  his  Eq.  (3).  Dougall fs  third  basic  solution  is  given  by 


u 

a 


eag^’3 


u  =  0 

Z 


*PZ -ease  note  that  lower  ease  Greek  letters  <t>3  0_,  ^  will  be  used  for  the 
Dougall  potentials  and  upper  case  ¥  for  the  Muki  potentials .  Such 
notation  is  intended  to  conform  with  that  in  References  1  and  and 
thus  makes  it  easy  for  the  reader  to  reconcile  the  results  here  with 
those  in  the  references .  However the  reader  is  advised  to  be  careful 
in  distinguishing  the  upper  and  lower  case  letters . 
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and 


which  are  his  Eqs.  (4).  As  explained  in  Appendix  B,  e  =  e 

n  otp  pa 

e12  =  X'  ■ 

In  Section  5  below,  we  shall  give  a  slightly  simplified  form  of  the 
Dougall  potentials. 

IV.  THE  MUKI  POTENTIALS 


r, 


Muki  expresses  his  potentials  in  terms  of  cylindrical  coordinates 
0,  z.  He  writes 


_  'Ll  ,  2  ay 

Ur  3r3z  r  30 


(11) 


u 


0 


I  3  $  2  — 

r  903z  “  3r 


and 


u 

z 


2(1  -  v)  V2$ 


(12) 

(13) 


where  $  is  an  arbitrary  biharmonic  function  and  ¥  is  an  arbitrary 
harmonic  function,  i.e., 


V 


n  =  o 


and  u  ,  u  ,  u  are  the  physical  components  of  displacement.  Our  Eqs. 
r  w  z  _ 

(11),  (12),  and  (13)  are  Mukifs  Eqs.  (3)  .  Writing 


$  = 


£ 

2A 


2£_ 

A 


and  using  Eq.  (6),  we  can  express  the  Muki  solution  in  cartesian 
coordinates  as  follows: 


u 

a 


2A  $ ,  + 

az 


eae  M’l 


u  =  2A  -  v2$ 

z  zz 


(14) 

(15) 


We  shall  find  this  slight  modification  of  the  form  of  the  equations 
to  be  useful  in  our  comparisons. 
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V.  HARMONIC  ANTI -DERIVATIVES  OF  HARMONIC  FUNCTIONS  - 
A  SIMPLIFIED  FORM  OF  THE  DOUGALL  POTENTIALS 

We  shall  demonstrate  here  two  closely  related  results  which  we 
shall  need  below: 

a.  If  g  is  a  function  such  that  g,^  is  harmonic,  then  there 
exists  a  harmonic  function  g  such  that  g,z  =  g,z- 

b.  If  f,  are  harmonic  functions,  a  =  1,2,  then  there  exists 

cx~ 

a  harmonic  function  f  such  that  f,  =  f,  ,  a  =  1,2 

’a  a 

We  shall  first  establish  b.  Given  that  f,  is  harmonic,  it  follows 
that  a 


tv  f),a  ■  o 


so  that  V  f  is  a  function  of  z  alone.  Let  h(z)  satisfy 


„ 2U,  ,  d2h 
V  h(z)  =  — j 

dz^" 


V2f 


then 

v2(f  -  h)  =  0  . 

(16) 

Let 

f  =  f  -  h 

(17) 

then 

it  follows  from  Eq.  (16) 

that  f  is  harmonic. 

Furthermore, 

since  h 

is  a 

function  of  z  alone,  h, 

a 

=  0,  whence  from  Eq. 

(17) 

which  establishes  b. 

The  proof  of  a  is  quite  similar,  but  we  shall  go  through  it  anyway. 

2  2 

Since  g,^  is  harmonic,  (V  g) ,  =0  and  thus  V  g  is  a  function  only 

of  and  x2 .  Let  £(Xj,  x2)  satisfy 

2  2 

V  %  =  i,  =  V  g 

aa  6  ’ 
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then 


V2  (g  -  Jl)  =  0 


(18) 


Define  g  by 


g  =  g  -  ^  ;  (19) 

then  by  Eq.  (18),  g  is  harmonic.  Moreover,  since  SL  does  not  depend  on  z 
i.e.,  £,  =  0,  it  follows  from  Eq.  (19)  that 

g  * z  =  g,z 


which  establishes  a. 

We  are  now  ready  for  a  modified  and  somewhat  simplified  form  of  the 
Dougall  potentials.  First,  let  <j>  be  given  by 


(20) 


Then  it  follows  that  <J>  is  harmonic.  Moreover,  given  any  harmonic 
<f>,  it  follows  from  the  immediately  preceding  argument  that  one  may  find 
a  harmonic  <j>  satisfying  Eq.  (20).  Thus  we  may  replace  Dougall!s  first 
solution  given  by  Eqs .  (10)  by 


u  =  Ax  6,  u ,  =  -  Ax  6,  -  d> 

y  yY,z  ’z  yY,y  y 


where  (j>  is  an  arbitrary  harmonic  function.  For  convenience,  we  shall 
write 


so  that  the  general  Dougall  solution  will  be  given  by 


u 

a 


Ax  6,  + 

aY '  z 


A0, 


*  Aeo.e*'e 


(21) 
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(22) 


Now  Eqs.  (21)  and  (22)  are  written  in  such  a  way  as  to  make  comparison 
with  the  Muki  solution,  Eqs.  (14)  and  (15),  easier. 

VI.  FINDING  THE  DOUGALL  POTENTIALS  GIVEN  THE  MUKI  POTENTIALS 

Suppose  that  a  solution  is  given  in  terms  of  the  Muki  potentials. 

We  address  ourselves  now  to  the  question  of  determining  a  corresponding 
set  of  Dougall  potentials. 

Because  the  second  term  in  Eqs.  (14)  has  the  same  form  as  the  third  term 
term  in  Eqs.  (21),  we  see  that  we  need  merely  find  the  Dougall  potentials 
<£,  0,  ty,  which  satisfy 

Axj,z  +  A0,a  +  Aea^,e  =  2A*,az  ,  (23) 


(24) 


given  a  biharmonic  Muki  potential  $.  Once  we  have  found  these,  we  may 
add  to  the  Dougall  solution  terms  of  the  form  Ae  to  get  agreement 
with  the  general  Muki  solution.  ^ 


If  we  calculate  A  from  the  Dougall  solution  with  the  use  of  Eqs. 
(21) ,  (22)  and  (9) ,  we  obtain 


A  =  (2A  -  1)  <f>,z 


2  2 

where  we  have  made  use  ofV<j>  =  V0  =  O.  Similarly  if  we  calculate  A 
from  the  Muki  solution,  as  given  by  Eqs.  (14)  and  (15),  we  obtain 


A  =  (2A  -  1)  (V2$), 


z 


4 

where  we  have  used  Eq.  (9)  and  V  4>  =  0.  Thus  4>  must  satisfy 


(25) 
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We  shall  see  that  it  will  suffice  to  take 


<f>  =  V  $  . 


(26) 


We  shall  discuss  below  how  the  Dougall  solution  will  be  altered  by 
replacing  <J>  as  given  by  Eq.  (26)  by  another  solution  of  Eq.  (25)  .  But 
first  let  us  pursue  the  implications  of  Eq.  (26). 

Substituting  Eq.  (26)  into  Eqs.  (23)  and  (24),  we  see  that  we  must 
now  find  0  and  ij>  such  that 


x  (V  4>) ,  +0,  +  e  \ p,  =  2$, 

a *■  J,z  a  afj  B  a  z 


(27) 


-x  (V  $) ,  +6,  =2$, 

or  J’a  *z  zz 


(28) 


We  start  with  Eq.  (28).  Let  0  be  any  solution  of  Eq.  (28)  for  0, 
i.e.,  let  0,^  satisfy 


0,  =  2$,  +  x  V  $, 

z  zz  a  a 


(29) 


We  shall  show;  then,  that  0,  is  harmonic.  To  this  end  we  differentiate 
Eq.  (29)  with  respect  to  x^.  (Remember  xa  g  =  6^)  .  We  obtain 


0,  „  =  2$,  „  +  6  vV  +  x  V2$,  „ 

zg  zzg  otg  a  a  ag 


2  2 

-  2$,  +  V  +  x  V  $ 

'zzg  g  a  'ag 


We  differentiate  again  with  respect  to  g  and  sum  on  g  to  get 


~  2  2  2 
0  =  2$  +V$  +6  +xV$ 

’zBB  ’zzBB  *BB  c*B  ’ctB  a  ’aBB 


2$,zzBB  +  2V  $,BB  +  XccV  $,aBB 


(30) 
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Now  we  differentiate  Eq.  (29)  twice  with  respect  to  z  to  get 


0,  =  2$,  +  x  V 

zzz  zzzz  a  azz 


Adding  Eq.  (30)  to  Eq.  (31),  we  arrive  at 


V  0,  =  0,  00  +  0, 

z  zgg  ’zzz 


=  2$  +  2$  +  2V  $ 

’zzgg  ’zzzz  *33 


+  x  [V2$,..  +  V2$,  ], 

a  L  SB  zzJ  a 


Note  that 


*,  oa  +  *,  =  V2$, 

’zzgg  ’zzzz  ’zz 


(v2$),zz  +  (v2$),gg  =  v2v2$  =  v4$ 


(31) 


(32) 


(33) 

(34) 


Putting  Eqs.  (33)  and  (34)  into  Eq.  (32),  we  get 

2"  4  4 

V  0,  =2  V  $  +  x  V 

z  a  a 


(35) 


2" 

But  $  is  biharmonic.  Thus  Eq.  (35)  gives  V  0,^  =  0. 

According  to  Section  5  above,  we  can  find  a  harmonic  function  0 

such  that  0,  =  0,  .  This  function  will  then  satisfy  Eq.  (28)  as  well 

2  z  z 

as  V  0  =  0.  It  will  be  our  second  Dougall  potential. 

We  have  only  left  now  to  find  a  harmonic  ^  which  satisfies  Eqs.  (27). 

We  shall  proceed  as  follows:  first,  let  ^  be  a  vector  (not  necessarily 

p 

the  derivative  of  a  scalar),  which  satisfies  Eqs.  (27).  We  shall  show 

that  ip  =  \\)  0.  This  will  then  allow  us  to  write  \pn  =  il)_  where  ij)  is  a 

.  p  > Y  p  ~  Bp 

scalar.  Then  we  shall  show  that  ip,  ^  is  harmonic.  Using  the  results  of 

Section  5  above,  we  can  then  assert  that  we  can  find  a  harmonic  function 
ip  which  satisfies  Eqs.  (27). 
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We  start  then,  with  ipa  satisfying  Eqs .  (27).  (A  cursory  examination 

P 

of  Eqs.  (27)  shows  that  and  \p2  are  determined  explicitly  once  $  and  0 

are  given.)  Then  let  us  differentiate  Eqs.  (27)  with  respect  to  x^  and 
sum  on  a.  We  get 


6  V2<I>, 

aa  z 


2 

x  V 

a  za 


aa 


e  ,  =  2$, 

a3  3  a  aaz 


or 


2V2<I>,  +  x  V2<I>,  +0,  +  e  ip  ,  =  2$,  . 

z  a  za  aa  a3  3  a  aa^ 


(36) 


Differentiation  of  Eq.  (28)  with  respect  to  z  gives 


-x  V  $  +  0,  =  2$, 

a  az  zz  zzz 


(37) 


Adding  Eq.  (36)  to  Eq.  (37),  we  get 


2V2$,  +  V20  +  e  „iK,  =  2V2$, 

’z  a8  p  a  ’z 


(38) 


Since,  however,  V  6  =  0,  Eq.  (38)  becomes  simply  e  ^  ,  =  0 

a  p  p  a 

or  ^2  1  "  ^1  2  =  ^  whence,  by  Green’s  theorem  in  the  plane,  there  exists 

a  scalar  \p  such  that  ip  =  ip,  ,  and  we  have  for  Eqs.  (27) 

a  a 


x  V  $ ,  +  0,  +  e  =  2$, 

a  *z  a  a3  3  az 


(39) 


Our  next  step  is  to  show  that  ip,0  is  harmonic.  We  differentiate 

3 

Eqs.  (39)  with  respect  to  x^  to  obtain 


6  V^O,  +  x  7^0,  +0,  +  e  nip,0  =  2$, 

ay  z  a  yz  ay  a3  3y  ayz 
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We  differentiate  once  more  with  respect  to  x^  and  sum  on  y.  We  get 

2  2 

25  V  0,  +  x  V  0,  +0.  +  e  =  20, 

ay  zy  a  zyy  ’ayy  aB  BYY  ’ayzy 


or 


2  2 

2V  0,  +  x  V  0,  +6,  +  e  iK„  =  20, 

az  a  zyy  ’ayy  afj  Byy  azyy 


(40) 


We  now  differentiate  Eqs,  (39)  twice  with  respect  to  z  to  get 


x  V  0,  +6,  +  e  J>,„  =  20. 

a  ’zzz  azz  aB  Bzz  ’azzz 


(41) 


We  add  Eqs.  (40)  to  Eqs.  (41)  to  obtain 


2V20 


.  +  x  V20,  .+  V20,  , 

az  a  I  yy  *zz 

=  2  fo,  +0,  *1  , 

yy  zz  (  ’a 


2  2~ 

+  V0,  +  e  V  jb.  _ 

a  aB  r,B 


Since  V  0,  =  0,  we  have 

ot 


2V20,  +  x  V40,  +  e  V2ti>,  =  2V20, 

az  a  *z  aB  B  az 


whence,  since  v  0  =  0, 


eaBV  *'B  =  °  ’ 


(42) 


,2? 


For  a  =  1,  Eqs.  (42)  give  V  =  0.  For  a  =  2,  Eqs.  (42)  gi 


give 


V  ip ,  ^  =  0.  Therefore 


V  =  0. 
a 


(43) 


15 


According  to  Section  5,  then,  Eqs .  (43)  imply  that  we  can  find  a 
harmonic  function  ip  such  that 


Such  a  function  will  be  our  third  Dougall  potential  and  we  are  done. 

Remark :  There  is  a  certain  amount  of  non-uniqueness  in  the  determination 
of  the  Dougall  potentials  from  the  Muki  potentials.  This  non-uniqueness 
we  shall  now  discuss. 

To  begin  with,  let  us  note  that  Eq.  (26)  is  not  the  general  solution 
of  Eq.  (25).  Indeed,  the  general  solution,  which  we  denote  by  is 

<t>  =  V2$  +  h(xa)  ,  (44) 


where  h(x  )  is  a  function  of  x, ,  x»  such  that 

CX  _L  Z 


V2h  =  h,  =  0 
aa 


(45) 


We  reserve  the  unbarred  quantities  <J> ,  0,  ip,  for  the  solutions  we 
already  have.  Suppose  that  0  and  \p  are  the  corresponding  other  Dougall 
potentials.  We  have  from  Eqs.  (23),  (24),  (44)  and  (45) 


2 

xV$.  +  0,  +  e  = 

a  ’z  ’a  agy,0 

2\z  • 

(46) 

-x  (V2$) ,  -  x  h,  +0, 

a v  J  ’a  a  a  z 

X  =  2$»zz  • 

(47) 

Comparing  Eqs.  (46)  and  (47)  with  Eqs.  (27)  and  (28),  we  get 


"  e)v+  =  0  >  (48) 

-xahJa  Me  -  e),z  =  0  •  (49) 
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It  follows  from  Eq.  (49)  that 


0=0+  zxgh,g  +  5(xp) 


Furthermore, 


V2z 


V'e 


2zh, 


YY 


zxQh,D 
3  BYY 


(50) 


It  follows  then  that  =  0.  We  put  this  result  into  Eqs .  (48)  to 

obtain 


(W’a  Z  *  +  W*  -  W>6  '  °' 


Now  both  xgh,g  and  £  are  plane  harmonic  functions.  Therefore  they  have 

harmonic  conjugates.  Let  n(x  )  be  the  harmonic  conjugate  of  xah,D,  i.e., 

a  p  p 


(x,  h,  ) ,  =  e  -ii,  . 

t  x  a  ag  6 


Let  £  be  the  harmonic  conjugate  of  then 

eag  tzn  +  C  +  *  "  ^>g  =  0  • 

It  follows  that 


i|j  =  i|j-zn-£;  +  az  +  b 


(51) 


where  a  and  b  are  arbitrary  constants.  This  is  true  because  ip  -  ip  +  zri  -  £ 
is  a  harmonic  function  of  z  alone.  Thus  Eqs.  (44),  (50)  and  (51)  give 
the  arbitrariness  in  the  Dougall  potentials. 
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VII.  FINDING  THE  MUKI  POTENTIALS  FROM  THE  DOUGALL  POTENTIALS 

Suppose  that  <p,  0  and  ip  are  a  set  of  Dougall  potentials.  Using 
Eqs .  (14),  (15),  (21)  and  (22),  we  see  that  a  corresponding  set  of  Muki 
potentials  $  and  V  must  satisfy 


Ax«+'z  *  M‘a  "  Aea6%  "  2M’cz  +  e«8A,,’6  • 


'AxJ’„  -  ♦  +  Ae>,  =  2A*,„  -  V  $ 

^  tx  Z  Z  Z 


(52) 

(53) 


Again,  we  equate  dilatations  to  obtain  Eq.  (25). 


Let  $°  be  a  particular  solution  of. 

V2$°  =  <(> 

Then  Eqs.  (25)  and  (26)  give 
V2($  -  *°),z  =  0  , 

2 

whence  V  ($  -  $°)  =  Ah(x  ),  h,  =0 

a.  ,aa 

Let  H,ctct  =  h-  Then  v2t$  -  -  AH]  =  0.  Thus, 

$  =  4°  +  AH  +  F (x. ) 


(54)' 


where  F(x  )  is  harmonic.  The  general  solution  of  Eq.  (25),  then  is  given 
by  Eq.  (54),  where  $  is  a  particular  solution  of  Eq.  (25),  F  is  an 
arbitrary  harmonic  function  and  H  is  an  arbitrary  biharmonic  function  of 
xa-  Let  us  Put  Ecls’  (26)  and  (54)  into  Eqs.  (52)  and  (53).  We  get 
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+  2F 


x  V2$°,  +  6,  +  e  „ip,„  =  2$>°,  +  2F,  +  e  V , 

a  .’z  'a  agr,g  ’oz  ’az  ag  'g 


-x  V2<J>°,  +  0,  =  2<J>°,  +  2F,  -  h 

a  ’ a  *z  *zz  *  zz 


or 


*  S"a  -  Zh,a  *  ea6  [*  '  T>’e  *  2*°-« 


-x  v2$°,  +  6,  =  2<J>°, 

a  ’a  *z  zz 


where 


(55) 

(56) 


6=0-  2F,z  +  zh 

(57) 

Now  let  n  be  the  conjugate  harmonic  function  to  h; 
Riemann  equations 

i . e .  ,  the  Cauchy 

h’a  =  eagn’g 

are  satisfied. 

Then  Eqs.  (55)  and  (56)  become 

x  V2$°,  +  6,  +  e  -i,.  =  2$°.  , 

a  ’z  ’a  agr’g  ’az  ’ 

(58) 

-x  V2$°,  +  0,  =  2$°, 

a  a  z  ’zz 

(59) 

where  if  =  \p  -  f  -  zn 

(60) 

Then  6  and  if  are  determined  as  in  Section  6. 
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Note  that  Eq.  (59)  determines  0  to  within  an  arbitrary  harmonic 
function  of  x^.  If  we  denote  this  function  by  P(xa)  and  a  particular 

solution  of  Eq.  (59)  by  0° 

0  =  0°  +  p(x  ) 
or 


If  q  is  the  conjugate  harmonic  function  of  p;  i.e.,  p>a  = 

and  if  ij;0  is  a  particular  solution  of  Eqs.  (58)  for  p  =  0,  we  have 

ij;  =  ijj0-q+az  +  b  , 


where  a  and  b  are  arbitrary  constants. 

We  now  go  back  to  Eqs.  (57)  and  (60).  We  have 

0°  -  0  =  -  2F,z  +  zh  -  p  (61) 

ip  -  =  t-  q-az-b  +  zn  (62) 

For  any  arbitrary  assignment  of  the  plane  harmonic  function  q,  and 
any  assignment  of  a  plane  biharmonic  function  H  and  any  constants  a  and 
b  we  can  determine  harmonic  functions  F,  and  ¥  from  Eqs.  (61)  and  (62). 
Then  F  is  determined  to  within  addition  of  an  arbitrary  plane  harmonic 
function  and  thus  $  is  determined  from  Eq.  (59) . 

It  follows  then  that  we  can  find  a  harmonic  function  ^  such  that 

=.  iJj,  and  thus  we  have  our  third  Dougall  potential.  In  this  way 
p  p 

we  see  that  one  can  always  find  harmonic  functions  <|>,  0,  which 
satisfy  Eqs.  (21)  and  (22)  whenever  u^  is  a  solution  of  Eqs.  (1). 
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APPENXIX  A 

NOTE  ON  CARTESIAN  TENSOR  NOTATION 


One  expresses  vectors  and  tensors  in  terms  of  their  components  in  a 
cartesian  coordinate  system  whose  coordinates  are  x^,  x^,  One  writes 

x^  for  the  general  coordinate  and  understands  that  i  will  range  over  the 

values  1,  2,  3. 

In  case  a  term  is  written  with  a  repeated  index,  summation  over 
that  index  is  understood.  Thus, 


x.v.  means  V*  x.v.  or  x,vn  +  x0v0  +  x_v7 
ii  11  22  33 

i  =  l 


a...  means  +  a.no  +'a.__ 

ijj  ill  i22  i33 


Generally  speaking,  the  index  on  a  vector  or  tensor  indicates  its 
cartesian  component.  Thus  v^,  v^,  v^  are  respectively  the  x^,  x^,  and 

x  components  of  a  vector  v. 

«J  A, 


A  comma  denotes  partial  differentiation  with  respect  to  one  of  the 
coordinates .  Thus , 


3v. 

i 


The  Kroneker  delta  represents  the  components  of  the  unit  matrix: 


6.  .  =  1  if  i  =  j  ;  6.  .  =  0  ifi^j  ;  and 

ij  ij 

we  have  6.  .  =  3. 

ii 

We  adopt  a  special  convention:  small  Greek  letters  range  over  the 
values  1,  2.  Thus,  u^  stands  for  u^  and  ^  only. 


xu  =  x.u.  +  x0u0, 
a  a  11  2  2 


Small  Latin  letters  continue  to  range  over  the  values  1,2,3. 
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The  quantity  6^g  is  the  two  dimensional  Kroneker  delta. 

represents  a  2  x  2  unit  matrix).  Thus,  6  =2. 

aa 

We  also  make  use  of  a  two-dimensional  alternating  tensor 
is  defined  by  the  following:  eag  =  -e^,  en  =  1.  Thus,  en 
eH  =  Similarly  e22  =  0.  Also,  e21  =  -1.  Its  matrix  is 

(: ;) 

To  see  the  use  of  this,  suppose  given  a  vector  v^  and  consider 

Wa  -  ea3v3’  ^us>  “  v2»  w2  =  ”vi'  This  is  used,  for  example,  in  the 

condition  that  a  two-dimensional  vector  be  the  gradient  of  a  scalar.  This 
condition  is 


3v. 

1 

2 

9x0 

3x, 

2 

1 

°r  Vl,2  = 

v9  , .  This  is  written  e  „v„  =  0 

^>1  a3  3, a 

Note 

that 

e  _e  =  6. 
a3  ay  3y 

This  can  be  checked  merely  by  trying  all  four  combinations  of  3,  y. 
For  example,  if  3  =  1,  y  =  2,  we  get 


(It 


e  .  It 
a3 


-e 


11 


Galea2 

’  ellei2 

+  e2ie22  =  ° 

> 

since  = 

e22  =  °* 

if  e  =  1,  y 

=  1,  we  get 

e  ie  i 

al  al 

=  ellell 

C12el2  =  ll 

since 

611  =  °* 

e12  =  lm 

Also  the  Cauchy-Rieman  equations  are  written  using  e  .  Thus  if  a 
function  f(z)  of  the  complex  variable  z  is  written  in  term^  of  its  real 
and  imaginary  parts,  f (z)=U(x1 ,x2)+iV(x1 ,x2) .  Then  the  Cauchy-Riemann 

equations  U,..  =  V,„  and  U,9  =  V}1  become  simply  U,  =  e  V,  . 

1  z  ^1  r  1  a  a3  53 
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APPENDIX  B 


DOUGALL  POTENTIALS  AS  GENERAL  SOLUTIONS  OF  THE  ELASTICITY  EQUATIONS 

We  shall  show  here  that  the  Dougall  potentials  give  the  general 
solution  of  the  elasticity  equations  as  represented  by  Eqs.  (7)  and  (8). 

To  verify  that  Eqs.  (21)  and  (22)  satisfy  Eqs.  (7)  and  (8),  as  long 
as  4>,  0  and  are  harmonic,  is  a  matter  of  straight-forward  substitution. 
We  shall  now  address  ourselves  to  the  more  involved  question,  namely 
that  of  showing  that  every  solution  is  of  the  form  given  by  Eqs.  (21) 
and  (22). 

Suppose  then,  that  u.  is  a  solution  to  Eqs.  (7)  and  (8),  or  simply 
Eqs.  (5).  Then,  1 


Differentiate  with  respect  to  and  sum,  to  get 
V2A  =  u. , . . .  =  0 


Thus  A  =  u. , .  is  harmonic.  Let 

i*i 


such  that  <f>  is  harmonic,  and  let  0  =  0  be  a  solution  of  Eq.  (22).  Then 
differentiating  Eq.  (22),  we  get 


2A*,ee  -  AV20 


z 


z 


But 


*’ee  ’  ■+ 


=  _ Z_ 

zz  2 A  -  1 


2AA, 

*  z 


2* 

-  AV  0 


z  2A  -  1 


z 
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But,  using  (8),  we  get 


2AA,z  2AA,z 

2A  -  1  =  2A  -  1 


2 A 

Av  0 , 


V  0,  =  0. 

*  Z 


Thus  we  can  find  a  harmonic  solution  0  of  Eq.  (22).  Now  let  \p  satisfy 

S 


u  =  Ax  <j>,  +A0,  +  Ae  \p 

a  a  z  a  ciB^B 


Differentiate  with  respect  to  a  and  sum  to  get 


u  ,  =  2A<j>,  +  Ax  <f> ,  +  A0,  +  Ae  n\p  , 

“  “  z  a  za  aa  ciB^B  a 


Let  us  then  differentiate  Eq.  (22)  with  respect  to  z  and  add  the  result 
to  this  equation.  We  get. 


(2A  l)*,z  +  AeaB*p,o  . 


Now  from  Eqs .  (B-l)  and  (B-2),  we  get 


A  =  A  +  Ae^iK, 

a  ' 


or 


(B-2) 


Thus,  =  i|),b  where  ip  is  a  scalar.  Thus  ip  satisfies  Eqs.  (21).  We  have 
just  to  show,  now,  that  ip , D  is  harmonic. 

p 
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From  Eqs.  (21),  we  get 


u  .  =  2A(f>,  +  Ax  <f> ,  +  A0,  +  Ae  ip, 

af  yy  za  a  zyy  ayy  &YY 


and 


u  ,  =  Ax  <}>,  +  A0,  +  Ae  01>,0_ 

a* zz  a  zzz  azz  ag  gzz 


Adding  the  last  two  equations,  we  have 


u  =  2A<f>,  +  Ae  V^i| i, 

a  T,za  ag  T’g 


Using  Eqs.  (7)  and  (B-l) ,  we  get 


2A  2 A  .  .  „2; 

A,  =  ™ - T  A>~  +  Ae  „V  if),, 


2A  -  1  a  2A  -  1  ’a  ag  r,g 


or 

2" 

V  *,g  =  0 

It  follows  then  that  we  can  find  a  harmonic  function  such  that 

^>o  =  ^ > o  and  thus  we  have  our  third  Dougall  potential, 
g  g 

In  this  way  we  see  that  one  can  always  find  harmonic  functions 
4> , 0 , ip  which  satisfy  Eqs.  (21)  and  (22)  whenever  u.  is  a  solution  of  Eqs. 
(1).  1 
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